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Simplifying products of fractional powers of powers 
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Abstract 

Most computer algebra systems incorrectly simplify 
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to rather than to 0/0. The reasons for this are: 

1. The default simplification doesn't succeed in simplifying the denominator to 
0. 

2. There is a rule that is the result of divided by anything that doesn't 
simplify to either or 0/0. 

Many of these systems have more powerful optional transformation and general 
purpose simplification functions. However that is unlikely to help this example even 
if one of those functions can simplify the denominator to 0, because the input to 
those functions is the result of default simplification, which has already incorrectly 
simplified the overall ratio to 0. Try it on your computer algebra systems! 

This article describes how to simplify products of the form w'^ (zi;^^) ^ • • • (w''") 
correctly and well, where w is any real or complex expression and the exponents 
are rational numbers. 

It might seem that correct good simplification of such a restrictive expression 
class must already be published and/or built into at least one widely used computer- 
algebra system, but apparently this issue has been overlooked. Default and relevant 
optional simplification was tested with 86 examples on 5 systems with n = 1. Using 
a spectrum from the most serious flaw being a result that is not equivalent to the 
input somewhere to the least serious being not rationalizing a denominator when 
that doesn't cause a more serious flaw, the overall percentage of most flaw types is 
alarming: 
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1 Introduction 

"When you are right you cannot be too radical" 

- Martin Luther King Jr. 

First, a few crucial definitions: 
Definition. Default simplification is what a computer-algebra system does to a stan- 



dard mathematical expression when the user presses ENTER or SHIFT ENTER , using 
factory-default mode settings without enclosing the expression in an optional transfor- 
mational function such as expand(. . .), factor(. . .), or simplify(. . .). 

Default simplification is the minimal set of transformations that a system does rou- 
tinely. Default simplification is called evaluation in Mathematica® and in some other 
systems. Any fixed set of default transformations is likely to omit ones that are wanted 
in some situations and to include ones that are unwanted in other situations. Therefore: 



• 



Most systems also provide optional transformations done by a function such as 
expand (...) or by assigning a certain value to a control variable such as trigExpand ^ 
true. 

• Some systems provide a way to disable default transformations. For example the 
Maxima assignment simp : false suppresses most simplification, whereas the Max- 
ima box(. . .), Mathematica Hold[. . .] and Maple freeze(. . .) functions suppress most 
or all transformations on their argument. 

Definition. Simplification is idempotent for a class of input expressions if simplification 
of the result (by the same default or optional transformations) yields the same result. 

Definition. A conveniently cancelable singularity is a removable singularity that 
can be removed exactly by functional identities such as sin(2u;) = 2 sin(tf;) cos(tf;) together 
with transformations such as a common denominator followed by factoring out the gcd of 
any resulting numerator and denominator, then using the law of exponents w^w'^ — )■ w^^'^ . 

For example, z'^z~'^ — )■ 2, sin(22;)/sin(z) — )■ 2cos(2;), and 
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which cancels the removable singularity at c = 0, leaving the non- removable singularity 
along the hyperbola ex = 1. However the removable singularity in sm[z)/z is not conve- 
niently cancelable because it can't be canceled exactly except inconveniently by means 
such as introducing the piecewise function 
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, otherwise. 



or the infinite series 



sin(z) >^ {—l)''z'^^ 
^^^(2^TI)!- 



Definition. A nested power product is an expression or a sub-expression of tlie form 

w^^iw^^y-'-iw^-y^ (1) 

witli n > 1, rational exponents, and a possibly or 1. 

This article describes simple algorithms that can be used in default and/or optional 
transformations to simplify nested power products correctly and well. The abstract 
presents one example of why this is important. 

Default and relevant optional transformations for Derive® 6.00, TI-CAS version S.lQj, 
Maxima 5.24.0, Maple*™ 15.00 and Mathematica 8.0.4.0 were tested on 86 examples for 
the simplest case where n = 1. The table in the Abstract shows the overall percentages 
for each of eight different decreasingly serious flaw types described in Section [3l 

Those large percentages for the six most serious kinds of flaws are alarming, and so 
are many corresponding percentages for each of the five systemso Wikipedia currently 
lists 29 other computer algebra systems, and I strongly suspect that most or all of them 
also have substantial room for improvement in this regard. 

Here is an outline of the rest of the article: Section [2] defines three more crucial terms. 
Section |3] describes eight prioritized goals for results that are nested power products, why 
they are important, and the reasons for their priorities. Section H] describes the tables of 
results at the end of this article and how the listed result flaws were measured. Section 
|5] describes four good forms for nested power products and how to obtain them: 

1. Form 1 merely standardizes the outer fractional exponents to the interval (—1, 1) in 
a way that doesn't introduce removable singularities, but instead tends to reduce 
their magnitude - perhaps completely. 

2. Form 2 further reduces many outer fractional exponents to [—1/2, 1/2] in a way that 
cancels as much of any removable singularity as can be done without resorting to 
form 4. Form 2 is an improvement on form 1 at the expense of more computation. 

3. Form 3 absorbs w" into one of the nested powers just prior to display if w" can 
thus be totally absorbed, giving a result with one less factor. Form 3 is an aesthetic 
improvement on form 2 at the expense of more computation. 

4. Form 4 completely cancels any cancelable singularity and nicely collapses all of 
the exponents into a single unnested exponent. However, this form often entails a 
complicated unit magnitude piecewise constant factor that is -1 raised to a com- 
plicated exponent. Unsophisticated users might be baffled by this factor, and even 
sophisticated users might abhor the mess. However, this form must be addressed 
because it can occur in input, it is valuable for some purposes, and some computer 
algebra systems generate this form for some inputs. 



^The computer algebra embedded in a succession of TI handheld calculators, Windows and Macintosh 
computers has no name independent of the product names, the most recent of which is TI-Nspire'™. 
^I am guilty as a coauthor of Derive and Tl-computer algebra. 



Section [6] suggests how to extend the algorithms to recognize syntactically different but 
equivalent instances of w in nested power products and how to extend the algorithms to 
some kinds of non-numeric exponents. Section [7] is an overall summary. The Appendix 
lists about one page of Mathematica rewrite rules that implement most of the third result 
form. Tables of results and their flaw numbers for the five systems and for the rewrite 
rules are at the end of the article. 

2 More key definitions 

In this article: 

• Unless stated otherwise, an indeterminate is a variable that has no assigned value, 
rather than a result such as 0/0. 

• Any finite or infinite-magnitude complex value can be substituted for indetermi- 
nates in expressions. 

• Fractional powers and square roots denote the principal branchjfl 

Definition. A canonical form for a class of expressions is one for which all equivalent 
expressions in the class are represented uniquely. 

Canonical forms help make cancellations of equivalent sub-expressions automatic. For ex- 
ample, if a computer algebra system always makes arguments of functional forms such as 
sin(. . .) canonical, then an input sub-expression such as sin ((x + 1)^) — sin [x^ -(- 2x -|- 1) 
automatically simplifies to rather than remaining unchanged as a bulky land mine that 
might make a subsequent result incorrect. Without canonical arguments, recognition 
of cryptically similar factors and terms requires costly tests such as determining if the 
difference in corresponding arguments can be simplified to 0. This might happen every 
time the same two functional forms meet during processes such as expansion of an in- 
teger power of a sum containing two sines, which can be often. In contrast, canonical 
arguments permit a much faster mere syntactic comparison of functional forms. 

As discussed in [H |6l [10] , canonical forms are unnecessarily costly and rigid for the 
entire class of expressions addressed by general-purpose computer algebra systems. How- 
ever, canonical forms are acceptable and good for default simplification of some simple 
classes of irrational sub-expressions such as nested power products. 

Definition. Zero-recognizing simplification for a class of expressions is simplification 
for which all expressions in the class equivalent to are transformed to 0. 

As illustrated by the example in the Abstract, a failure to recognize that a sub- 
expression is equivalent to can lead to dramatically incorrect results. Therefore it is 
desirable for default simplification to have at least a zero-recognition property. It has 



^By default some systems assume that indeterminates represent real values and/or use the real branch 
wherein for reduced integers m and n, (—1)™/" — > 1 for m even, and (—1)'"/" — > —1 for m and n odd. 
However, most computer algebra systems provide a way to force the principal branch if it isn't the 
default - and to declare that an indeterminate is complex if that isn't the default. 



been proven impossible to guarantee this even for some rather simple classes of irrational 
expressions, but a strong effort should be made to achieve at least zero recognition for 
as broad a class of expressions as is practical. 

Definition. Candid simplification produces results that are not equivalent to an ex- 
pression that visibly manifests a simpler expression class. 

For example, in a candid result there are no superfluous variables, degree magnitudes 
are not larger than necessary, there are no unnecessary irrational sub-expressions, and 
irrationalities are nested no more deeply than necessary. Thus without being as rigidly 
constrained as canonical forms, candid simplification yields more desirable properties 
than mere zero-recognizing simplification. 

Definition. In this article undefined means an unknown point in the entire infinite 
complex plane, such as the result of 0/0|j 

Definition. A conveniently representable subset of the infinite complex plane is one 
that is reasonably representable using constant expressions extended by sets, intervals 
and the symbol oo. 

Conveniently representable proper subsets of the infinite complex plane are regarded 
here as defined. Particular computer algebra systems might not be able to represent the 
full range of possibilities, but this article is suggesting what should be done as well as 
reporting the current situation. These ideas are discussed in more detail in [9], but for this 
article the major defined subset of interest that isn't a single point is the result of u/0 for 
any particular non-zero complex constant u. This result should be some representation 
of complex infinity. Among many other benefits it permits the correct computation 
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Does your computer algebra system do this? 

• For Mathematica, 1/0 —)■ Complexinf inity. 

• For Derive, with its default real domain, 1/0 —t- ±oo. 

• For TI-CAS, regrettably 1/0 -^ undef. 



"'it is of course audacious to define undefined. Although unnecessary for this article, systems could 
usefully also 

• display 0/0 as 0/0 rather than a vague controversial word such as "undefined", and 

• contract functions of 0/0 to strict subsets of the complex plane wherever possible, such as 
arg(0/0) -^ (— tTjTt]. Having arg(0/0) — >■ 0/0 snatches defeat from the jaws of compromise. 
Try this on your systems! Many systems throw an error, which is worse because it requires even 
amateur authors of functions to know about all the potential throws, catch them or vet to prevent 
them, and respond appropriately to make their functions robust. 



• Maxima and Maple inconveniently throw an error. 

When a proper subset of the infinite complex plane isn't conveniently representable, then 
the next best thing is to degrade it to 0/0. However, that shouldn't be done for subsets 
that are as easily represented as complex infinity. 

If finite or infinite magnitude complex numbers are substituted for all of the indeter- 
minates in an unsimplified input expression, then that input expression is undefined at 
that point if and only if the result is 0/0. 

Definition. A generalized limit is the set of uni-directional limits of an input expres- 
sion from all possible directions in the complex plane. 

When the generalized limit of an input expression at a conveniently cancelable singu- 
larity is a conveniently representable proper subset of the entire infinite complex plane, 
then this article regards it as not only acceptable but commendable to cancel the singu- 
larity and thereby produce a result expression whose substitutional value is that conve- 
niently representable subset at that point. 

Reasons for this attitude about mathematics software include: 

• Otherwise the results tend to be unacceptably complicated^ 

• There is a high likelihood that the physical problem is actually continuous there 
too - Nature abhors a removable singularity. Removable singularities are often an 
artifact of the modeling such as using a polar or spherical coordinate system. 

• Cancelable singularities are often a result of an unnecessary previous transfor- 
mation unavoidably done by a system (such as inappropriate rationalization of a 
denominator) or a result of a previous transformation such as monic normaliza- 
tion, a tangent half angle substitution, or expansion into partial fractions deemed 
necessary to obtain an anti-derivative. 

• Cancellation to simplify nested power products is consistent with quiet transforma- 
tions such as w/w — ;■ 1 that are currently unavoidable in most computer algebra 
systems; 

• Symbolic cancelation tends to reduce rounding errors near removable singularities 
for subsequent substitution of fioating-point numbers. 

However, this transformation of expressions has the composability consequence that sub- 
stitution of numeric values doesn't necessarily commute with simplification. To accom- 
modate either treatment of expressions, computer algebra systems could and should build 
in provisos such as "| u; ^ 0" that are optionally attached automatically to intermedi- 
ate and final results containing canceled removable singularities, as suggested in [3l [TO] . 
Meanwhile, implementers who do not want to completely cancel cancelable singularities 
for simplifying nested power products can adapt the algorithms presented here to merely 
reduce the magnitude of cancelable singularities, such as {z"^) j z^ — > {z^) / z rather 
than transforming all the way to z\fz^ . 



^Canceling a gcd occasionally increases bulk significantly, such as (a;^^ — l)/(a; — 1) ^x^^ + a:;^^ + --- + 
cc + 1, but the algorithms described here consider only syntactic cancellation, which always decreases 
bulk. 



3 A list of goals for simplifying nested power products 

The most important concern is correctness, followed by candidness, then aesthetics and 
compliance with custom. More specifically, here is a list of desirable but partially con- 
flicting goals for simplifying nested power product and their differences, in decreasing 
order of importance: 

1. The result should be equivalent to the input wherever the input is dehned. (It is 
acceptable for the result to be a generalized limit of the input where the input is 
0/0.) 

2. A linear combination of two or more equivalent nested power products should sim- 
plify to a multiple of a single power product - or to if the linear combination is 
equivalent to 0. 

3. Let the net exponent of (^w'^'^Y'' be 

Afc := /3fc7fc, 

and for a product of nested powers of w let the total positive nested exponent 
and the total negative nested exponent be 

n 

A+ := ^max(Afc,0), (2) 

fc=i 

n 

A_ := ^min(Afc,0). (3) 

fc=i 

When possible, use the transformation 

(w^^y ^ w'^'^f'' {w^^y'-""- (4) 

with appropriate integers nik to minimize min (max {a, 0) + A+, — min (a, 0) — A_), 
thus canceling as much of any removable singularity as is possible by this means. 

4. When possible, fuUy absorb the w°' into the nested powers of (w^''^ '' to have fewer 
factors. 

5. Otherwise use transformation (jl]) to minimize A+ — A_ to minimize the contribu- 
tions of the troublesome nested powers. 

6. Inputs that are equivalent where both are defined should produce the same (canon- 
ical) result. 

7. Results should be idempotent: Reapplying the same default or optional simplifica- 
tion to the result should leave it unchanged. 

8. To help achieve goal 6, rationalize a denominator in a nested power product when 
this doesn't introduce a removable singularity or increase its magnitude. 



A larger numbered goal should not be fulfilled if the only way to fulfill it is to violate a 
smaller numbered goal. For example, fulfillment of goals 5 or 8 can often violate goals 1, 
2 and/or 3. 

The reasons for this ranking of the goals are: 

1. A violation of goal 1 is most unsatisfactory because it is a result that is not equiv- 
alent to the input everywhere the input is defined. For example, if expression w 
can be 0, then rationalizing the denominator of l/y/w to give ^Jwjw makes an 
input that is a well-defined complex infinity at u; = become 0/0 there. A more 
serious example is the mal-transformation (t(;~^)~^/^ — )■ (u;^)^/^, because the two 
sides differ along the entire positive and negative imaginary axis. For example 
(i~^)~^/^ = — z, whereas ii^)^!'^ = i. 

2. The example in the Abstract shows the importance of zero recognition. For example 
if default simplification of one nested power product produces yfzj z and default 
simplification of another nested power product produces the equivalent expression 
z/ y/z, then the latter violates goal 8 and together they violate goal 6. These 
violations are minor; but if default simplification doesn't simplify their difference 
to 0, then that is a violation of goal 2, which is serious. 

3. A violation of goal 3 is next most serious because it is a squandered opportunity to 
improve the result by canceling a conveniently cancelable singularity and thereby 
making the result have the limiting value at w; = rather than be undefined there. 
For example, 

/ 2\5/2 2 / 2\3/2 

- — - — -^ ^ — - — -^ w (w^) , (5) 

w w ^ ' 



4. A violation of goal 4 is more complicated than need be. For example, most people 
would J 
factor. 



would agree that w^ (w^) is more complicated than (w^) , which has one less 



5. Goal 5 is important because when there is more than one factor of the form {w^^ '^ , 
there might be more than one way to distribute only some of w°' into the nested 
powers. In contrast if 7^ is a half-integer power then there are only two ways to 
minimize |7fc| by factoring an integer power of w" out of {w°')'^, or only one way 
for other fractional powers. Moreover, unnested exponents are less specialized and 
can therefore interact more freely with other factors in a product. For example, for 
intermediate results (jS]) and (^, 

^(^2)3/2 ^ ^2(^2)1/2^ 

1 / 1 y/' 1 / 1 ^'/' 

W \W'^ I U;2 \ y;2 



6. For a given w, the above goals tend to yield the most concise possible nested power 
product in terms of w. Therefore it is better to have the consistency of having 
two inputs that are equivalent where they are defined return the same most concise 
form. More importantly a canonical form for nested power product sub-expressions 
greatly facilitates achieving goal 2 - a major benefit for very little effort. 

7. Without idempotency, an unaware user could obtain inconsistent results, and a 
cautious aware user would have to re-enter such results as inputs until they cycle 
or stop changing. 

8. Goal 8 complies with the custom of rationalizing denominators and helps achieve 
canonicality goal 6. For example, 



■U;2 W'^ W 

But rationalization should not be done at the expense of lower-numbered goals. 
For example, 

1 ^ iw^f' 



^(^2)2/3 ^3 ' 

because although it reduces the absolute value of the outer exponent (goal 5), it 
violates goal 1 by making an input that is complex infinity at u; = become a 
result that is 0/0 there. 



4 Important information about the tables 

Tables at the end of this article show the results that occurred for each example with 
each system and with the Appendix rewrite rules. In all of the tables the goal numbers in 
the Section [3] goals that aren't satisfied but could be satisfied without violating a lower- 
numbered goal are listed beside each result. Unmet goal numbers 1 and 2 are boldface 
to emphasize their extreme seriousness. 

4.1 Examples, test protocol and table interpretation 

Tables [2] through [7] report default and relevant optional transformation results for test 
family 1: multiplying w;™ by (tf^)" for successive integer m = —3 through 3 in 
combination with successive integer n = —3 through 2. For comparison with results that 
meet all of the goals. Table [H has corresponding results for form 3 described in Section 
El as produced by the one page of Matheinatica rewrite rules listed in the Appendix. 

Tables [9] through [15] report default and relevant optional transformation results for 
test family 2: multiplying w™ by (w~^)" for the same combinations of m and n. 
For comparison with results that meet all of the goals. Table [H] has corresponding form 3 
results produced by the Appendix rewrite rules. 



To help assess compliance with goal 8, Table [T6] compares results for all of the systems 
with the Appendix rewrite rules on the particularly simple input w /yw^ . Some of the 
examples in test families 1 and 2 also test this goal. 

Table [T7] tests only whether or not the expression 



w 



-(-1) 



2farg(™ )-2arg(iu) j/n 



simplifies to 0. This is the difference between equivalent expressions in form 3 and form 
4. This is a more difficult but not impossible problem. All systems fail - including the 
Appendix rewrite rules, which do not address this issue. 
Here is how compliance with the goals was assessed: 

1. Most of the results that violated goal 1 did so only at w = 0. However, the goal 1 
violations in Tables [T3] and [T^ instead or also are not equivalent to the input where 
it is defined along the entire positive and negative imaginary axis. This is caused 
by an outlaw of exponents: transforming [w~^^ to [w^^ for fractional /i, which 
is not valid along these semi-axes. 

2. For test family 1, each input is equivalent to the input two rows down and one 
column left wherever both are defined, and their omnidirectional limits are iden- 
tical wherever one of the inputs is 0/0. Therefore to assess compliance with 
goal 2 (zero-recognition), for every entry in the table I computed the difference 
w'" {w'^)^ — iLi'"+2 {w'^Y o^ ^he optional transformation thereof and the 
difference w™ (w^)" — w™^^ {'^'^Y 5 then considered it a flaw for the entry 
if either of these two differences was non-zero. Thus compliance with this goal is 
not discernible from merely inspecting the result entries. Compliance is a property 
of the default simpliflcation or optional transformation when given the difference 
of two non-identical but equivalent nested power products. For test family 2, each 
input is equivalent to the input two rows down and one column right wherever 
both are deflned, so I did an analogous test for that. It is of course possible for 
an entry to pass these limited tests but fail for more widely separated equivalent 
inputsjj For Table [16] the result is equivalent to -Jw/w and w^/(w^)^/^, so I tested 
whether or not the corresponding two differences or optional transformation thereof 
simplifled to 0. Table [T71 tests zero recognition directly, and only than - but with 
only one particular difference in equivalent forms rather than only two. 



3. To comply with goal 3 without violating lower-numbered goals, a result w°' {w^^ 



7 



^This happens for Tables [13] and [14] All of the columns would exhibit flaw 2 if one of the two 
equivalent expressions was always taken from the correct results in columns 3 or 4. The results are not 
equivalent to the inputs for columns 1, 2, 5 and 6, so the only reason the difference simplified to for 
columns 1 and 6 was the subtraction of incorrect but identical results - an instance where two wrongs 
make a right. 
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should have 

a = V 

(sign (d) = sign (/37) A I7I < 1) V 

(a 7^ A sign (a) ^ sign (/37) A I7I < 1 A 

min(|/37| , |d|) < min(|/3 (7 - sign (7))! , |d + /3sign (7)!)). 

4. To comply with goal 4, a should not be an integer multiple of /3. 

5. To comply with goal 5 without violating lower numbered goals, 

d = V (sign(d)=sign(/37) A |7|<1) V I a 7^0 A sign(d)7^sign(/37) A |7|<- 

6. For compliance with goal 6, every result entry for test family 1 should be identical 
to the entry 2 rows down and one column left, whereas every result entry for test 
family 2 should be identical to the entry 2 rows down and one column right. To 
equally assess the top two rows, the bottom two rows, the leftmost column and 
the rightmost column, I computed extra neighbors bordering those shown. When 
there were differences, I did not penalize the best displayed results for equivalent 
entries unless their was a better displayed result one column and one or two rows 
outside the table. However, I did penalize all of the not-best members for equivalent 
entries. I similarly tested the result in Table [16] against the equivalent expressions 
y/w/w and w"^ / {w'^Y^'^ . It is of course possible for an entry to pass these tests but 
fail for more widely separated equivalent inputs. 

7. To test compliance with goal 7, I resimplified each result with either default sim- 
plification or the optional transformation used for the original input, then checked 
for identical results. Compliance with this goal is not discernible from merely in- 
specting the result entries. 

8. To test compliance with goal 8, I manually rationalized the results having a frac- 
tional power in the denominator and a 7^ by multiplying the numerator and 
denominator by (ty"^)^/^ for test family 1 or ynp- for test family 2. It was counted 
as a flaw if and only if that forced rationalization did not introduce a removable 
singularity or increase its magnitude. 

4.2 Remarks about particular results. 

Maxima also has a relevant rat(. . .) function. For these examples, it generally produces 
the same result as default simplification, except that fractional powers are represented 
as an integer power of a reciprocal power - or an integer power of JT7'. for half-integer 

powers. Thus a default result w^ {w"^) would instead be w'^{{w'^Y'^Y, and a default 
result w^ (w"^) would instead be w^yuP . The standard definition of m™/" for reduced 
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integers m and n is (m^/") , which is consistent with the alternate definition ein(M)m/n|7| 
Consequently, the Maxima rat (...) function makes the standard interpretation of the 
result more explicit at the expense of clutter. Nonetheless, it might be helpful as a 
precursor to semantically substituting a new expression for (w"^) via syntactic substi- 
tution in a expression containing (w"^)"^ for several different integer m. For the sake 
of brevity, results are not included for the rat(. . .) function because its flaws were very 
nearly identical to default simpliflcation, regarding ((w^)^/")'" as (w^)™/". 

Mathematica, Maxima and Maple also respectively have relevant Power Expand [...], 
radcan(...) and simplify(..., symbolic) functions. However, they always transform (w^)"^ 
to w^'^, which is not equivalent along entire rays from w = 0. I didn't test these functions 
because their purpose is presumably partly to allow these risky unconditional transactions 
for consenting adults. However, these three systems. Derive and TI computer algebra 
also have safe ways to enable such desired transformations, when justifled, by declaring, 
for example, that certain variables are real or positive. 

"... a man who thought he could somehow pull up the root without affecting the 

power." 
-adapted from Gilbert K. Chesterton 

To make sure that w is regarded as a complex variable and the principal branch is 
used rather than the real branch: 

• All of the Derive results follow a prior declaration w :G Complex. 

• All of the TI-CAS results necessarily used w_ rather than w to manifestly declare 
it as a complex indeterminate. However for consistency w_ is displayed in all of 
the tables as w because Table [12] is shared with Maple for brevity. 

• All of the Maxima results followed a prior assignment domain : complex and a prior 
declaration declare (w, complex). 

As illustrated by Tables [T^ through IT7\ the Maple simplify(. . .) function expresses half- 
integer powers of squares or of reciprocals of squares using the Maple csgn(. . .) function 
deflned by 

fl, if 3fJH > V 5R(w;) = A S^H > 0, 
csgn(ti;) := <^ ^ ' ^ ' v / - ' (^^^ 

1—1, otherwise. 

The right side of this deflnition is a simplifled special instance of form 4, for which csgn(u;) 
is a convenient abbreviation for those familiar with itjj 
Regarding Table [T71 



^This is not generally equivalent to (u™) '": "be faithful to your roots" - Mason Cooley. 

^Jeffrey [S] uses the unwinding function to generalize csgn to a C„ that works for all fractional powers. 
If and when implemented in Maple, that will avoid unwelcome mixtures of csgn(. . .) with other form 4 
notations for results containing both half-integer and other nested powers. 
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Mathematica (hence also the Appendix rewrite rules) did the automatic transfor- 
mation 

(_l)(V2)(Argp])-2ArgH)/7r _^ ^(Argp])-2ArgH)/7r^ 

Although it eliminates the 1/2 factor from the exponent in this case, it does so at 
the expense of candidness by introducing i into an expression that is real for all w. 

For TI-CAS the arg function is spelled "angle" and regrettably angle(O) is returned 
unchanged rather than transforming to 0. Therefore the input was 

V'^^_ _ /_-, Nwhen(ui_=0, 0, (l/2)(anglc(«)_2)-2angle(io_))/7r). 
W_ 

As indicated in Table [T71 the real power of -1 in the input was changed to an imag- 
inary power of e in the result. This has the candidness disadvantage of introducing 
i into an expression that is real for all real w_. 

For Derive the arg function is spelled "phase" and regrettably phase(O) returns 
7r/2 ± 7r/2, which denotes an unknown element of {0,7r}, which are the only two 
possibilities for real arguments. Therefore the input was 

""^^ _ / -|\IF(to=0, 0, (l/2)(PHASE(«;2)-2PHASE(io))/7r). 



W 

• For Maple, the arg function is spelled "argument", and for Maxima it is spelled 

"carg". 

If you are interested in results for some other systems, then try a few of the examples 
that are heavily flawed for most of the five tested systemso First do whatever is neces- 
sary so that fractional powers use the principal branch and w is regarded as a complex 
indeterminate. 



5 Four alternative forms 

Section [1] explains the reasons for four separate forms. 

5.1 Form 1: Reduction of outer fractional exponents to (-1, 1) 
Definition. For x G M the integer part function 

I t^\ .= ; W ' if ^ - 0' 
\x~\ , otherwise. 



^If you are familiar enough with those systems, then most of them probably have a quick way 
to generate all of the results for test families 1 and 2 by entries analogous to the following one for 
Mathematica: 

Table [Table K(u;2)'=,{fc, -7/3, 8/3}, {j, -3, 3} //TableForm 

I am interested in knowing your results. 
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Definition. For x G M the fractional part function Fp(x) := x — Ip(x). 
Proposition 1. For /9 G Q, 7 G Q — Z, and arbitrary expression w G C, 

(^/3)^^^/5Ip(7)(^/3)FpW. (8) 

Proof. We have 

(^/3)^^(^/3)IpW(^/3)FpW (9) 

because: 

1. With 7 G Q - Z, Ip(7) = V sign (Ip(7)) = sign (Fp(7)). 

2. For any expression m G C and ri,r2 G Q | ri = V signri = signr2, 

^ri+r2 ^^n^r2^ (10) 

even at m = with ri and r2 both negative, making both sides of (ITOj) be complex infinity. 

3. By Proposition m we also have (w^) = i^^^pW because Ip (7) G Z. D 

Therefore Form 1 is simply to transform w" (w;^^) ^ ■ ■ ■ (u;^") " toward canonicality 
by transforming every positive fraction 7^ to the interval (0, 1) and every negative fraction 
7fc to the interval (—1, 0). The various y7^p('>fc)^'= are combined with the original w°', giving 
a transformed expression 

W := w^w^^y' ■ ■ ■ {w^-y- (11) 

where a might be 0. 

This form 1 satisfies goals 1 and 7, while possibly contributing progress toward goals 
2, 3, 5 and 6. This form also has the advantage that if (i/;^) is subsequently raised to 

any power A, then we can simplify it to the simplified value of (w^) by Proposition [3] 
because — 1 < 7 < 1. For example, 



{w^f'f" ^ {w'f 



Although there is no such thing as a free radical in computer algebra, this transformation 
of each nested power is fast and easy to implement because it occurs only for certain 
fractional powers of powers, which are relatively rare, and very little work is done even 
when it does occur. There is no good reason why default simplification shouldn't do at 
least this much. 

However, default simplification for many systems unavoidably collects similar factors, 
resulting in a partial reversal of this transformation whenever a resulting unnested ex- 
ponent a is identical to one of the inner nested exponents. This happens for Derive, 
TI-CAS and Matheniatica, but not for Maple or Maxima. With unavoidable collection, 
unconditional magnitude reduction of fractional outer exponents can lead to an infinite 
recursion such as 

w (w^) -> w"^ (w^) -^ w [w'^y -)■■■•. 

Therefore in the Appendix Matheniatica rewrite rules: 
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• Transformation (w^) — )■ w^'^^"'^^ (w^) is used unconditionally only prior to 
default simplification. 

• The rewrite rules that are active during default and optional transformations do not 
reduce the magnitude of 7 if doing so would give an unnested exponent ^identical 

to/3. 

• This transformation is used after default simplification only if there is an unnested 
factor w"" and the transformation would not be reversed by unavoidable collection 
of similar powers. 

Implementations for other systems might have to overcome this difficulty in some other 
way, or compromise and not always produce a form with outer fractional exponents in 
the interval (—1, 1) when w°' can't be fully absorbed into some nested power. 

5.2 Form 2: Further reducing some outer exponents to (-1/2, 
1/2| 

Form 2 is form 1 supplemented by an additional transformation. 

Expression W given by definition fITT]) is equivalent to expression W everywhere that 
W is defined, because at the only questionable point w = 0: 

1. Expressions W and W are both if a > and all of the (!^klk are positive. 

2. Otherwise expression W and W are both complex infinity if a < and all of the 
/3fc7fc are negative. 

3. Otherwise if d > and all Pk^k > 0, then W is 0/0 but W has improved to 0. 

4. Otherwise if d < and all Pk'jk < 0, then W is 0/0 but W has improved to complex 
infinity. 

5. Otherwise both W and W are 0/0. However, the magnitude of the multiplicity of 
the removable singularity is less for W if for any 7^,-, |7a;| > 1- 

Expression W is canonical in cases 1 through 4, but not necessarily for case 5. For 
example, 

1. The different equivalent expressions 2~^ (2;^) and z (z"^) both have outer expo- 
nents in (-1, 1). Of these two alternatives, the latter is preferable for most purposes 
because the |— 2/3| < |4/3|, making multiplicity of the uncanceled portion of the 
removable singularity have a smaller magnitude. Thus a rationalized numerator is 
sometimes preferable to a rationalized denominator. 

2. The different expressions z~^ (z^) and z~^ (z^) both have outer exponents in 
(-1, 1), and they are equivalent wherever the first alternative is defined. However, 
the latter unrationalized denominator is preferable because the former is 0/0 at 
z = where the latter is defined and equal to the complex infinity limit of the 
former. 
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3. The different equivalent expressions z (z^)~ and z"^ (z^) both have outer expo- 
nents in (-1, 1), and the niultiphcities of the uncanceled portion of their removable 
singularity at 2; = are both 1. Of these two alternatives, the latter is slightly 
preferable because it has a traditionally rationalized denominator rather than a 
rationalized numerator. 

Thus after producing form 1 we can sometimes add 1 to a negative 7^ or subtract 1 from 
a positive 7^, then adjust a accordingly to reduce the magnitude of the overall removable 
singularity - perhaps entirely. If not, perhaps we can at least contribute toward goals 2, 
6 and 8 by rationalizing a square root in the denominator. 
Let 

Afe := /3fc7fc, 
A := a + Ai + --- + A„. 

Transforming any of the {w^^^ * to w"^^^^ {w^^^ ^ * for any integer m^ leaves A 
unchanged. 

Our primary goal is, whenever possible, to make all of the A/^. have the same sign and 
for a to have either the same sign or be 0. A secondary goal is to prefer —1/2 < 7^ < 1/2. 
Therefore, the algorithm to convert form 1 to form 2 is: 

1. If A > 0, then for each A^ < 0, add sign (/3fc) to 7^ and subtract |/3fe| from a, then 
return the result. 

2. If A < 0, then for each A^ > 0, subtract sign (/3fc) from 7^ and add |/3fc| to a, then 
return the result. 

3. For each 7^ > 1/2, subtract 1 from 7^ and add /3fc to a. 

4. For each 7 < —1/2, add 1 to 7^ and subtract /3fc from a. 

5. Return the result. 

This canonical form 2 satisfies all of the goals except for the aesthetic goal 4. 

For brevity the Appendix rewrite rules consider only one A^ at a time. This is 
sufficient for all of the test cases, which have only one nested power. For an industrial- 
strength implementation, each time we multiply a fractional power of a power by a 
product of one or more factors, we should inspect those factors for identical expressions 
w and apply the above algorithm if that subset if non-empty. The cost is O {ric) where Uc 
is the number of cofactors. The opportunity occurs only when multiplying a fractional 
power of a power, which is rare; and the number of factors in a product is typically quite 
small. Therefore it also quite reasonable to do this in default simplification. 

5.3 Form 3: Finally, fully absorb w" into a fractional power if 
possible 

Form 3 is form 2 followed by an additional transformation. 
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Form 2 can result in an expression such as z'^ (2;^) , for which many users would 
regard (2;^) as a simpler result because it has one less factor. We can often absorb at 
least some of 2° into one of the [z^^-^ ^ by the transformation 

which doesn't change the domain of definition. However, this transformation seems in- 
advisable unless Fp [a/ fik) = 0, because otherwise it increases the contribution of a 
troublesome nested power without reducing the number of factors. Also, this transfor- 
mation is problematic during intermediate computations even if Fp (a/Pk) = 0, because 
when there is more than one nested power, then more than one might be eligible, making 
it awkward to maintain canonicality achieved by form 2. Moreover, absorption conflicts 
with transformations done to obtain form 1 or 2, thus risking infinite recursion. 

A solution to this dilemma is to fully absorb w°' only just before display - after all 
other default and optional simplification. This does have the minor disadvantage that 
what the user sees doesn't faithfully represent the internal representation. However, that 
bridge has already been crossed by most systems, which for speed and implementation 
simplicity internally use, for example, (. . .)^/^ to represent a displayed ^/TTT and a-\ — 1*6 
to represent a displayed a — b. 

When there is more than one nested power of w, then there might be more than one 
way to absorb a completely into those nested powers. For example, 

6 / 2\l/2 / 3\l/2 / 4\l/2 / 2\V2 / 3\ 1/2 / A\l/2 

In general, the possible resulting expressions are given by 

where the tuple of integers {7711,7712, . . . , m„) is a solution to the linear Diophantine equa- 
tion 

mi/3i + 7712/92 H h 77lnl3n = Oi- 

Solutions exist if and only if a is an integer multiple of gcd (/3i, /32, . . . /3„), in which 
case there might be a countably infinite number of tuples. However, to avoid introducing 
removable singularities or increasing the magnitude of their multiplicity, we are only 
interested in solutions for which sign(mj/3j) = sign (a) for j = 1,2, ... ,n. Papp and 
Vizvari [7] describe an algorithm for solving such sign-constrained linear Diaphantine 
equations, and the Mathematica Reduce [. . .] function can solve such equations. For 
example, suppose our canonical form 2 result is 



z 



^^ (z'^/'rf {z^'/'f. (13) 
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In Mathematica, we can determine the family of integers rrii > to add to 1/2 and 
7712 > to add to 1/3 that together absorb z^'^ as follows: 



-77ii + —77i2 == 14 && -mi > && ^^^12 > 0, {nii,m2}, Integers 



ln[l] : = Reduce 

//TraditionalForm 

Out[l]//TraditionalForm = (?7ii = 3 A 7712 = 8) V (?tii = 8 A ?7i2 = 5) V (?7ii = 13 A ?7i2 = 2) 

Regarding the choice between alternative absorptions, canonicality is not as important 
for a final displayed result as it is during intermediate calculations where it facilitates 
important cancellations. However, with more than one solution, we could choose one in 
a canonical way as follows: Order the Pj in some canonical way, such as the way they 
order in (w^^^ ^ ■ ■ ■ (w^") ", then to choose the solution for which 7711 is smallest, with 
ties broken according to which ?7i2 is smallest, etc. 

Solution of sign-constrained linear Diophantine equations can be costly - probably 
too costly for default simplification. Consequently, the rewrite rules in the Appendix 
simply absorb if" if and only if it can be completely absorbed into a single power of 
a power, in which case the particular one is the first one encountered by the pattern 
matcher. This is canonical, but it doesn't absorb w°' for examples such as flT3|) . However, 
this transformation is inexpensive because it is done only once in one pass over the 
expression just prior to display, and the transformation requires comparing a power of a 
power with its cofactors only in products where powers of powers occur. 

All but this absorption rule are automatically applied before default simplification so 

that, for example, the input 

w — w 



correctly simplifies to indeterminate, meaning 0/0, rather than to 0. 

The rewrite rules in the Appendix are not much more than the minimal amount 
necessary to generate the form 3 results in Tables [1] and [HI together with the relevant 
rows in Tables [16] and [T71 

5.4 Form 4: One unnested power times a unit-magnitude factor 

Form 4 is quite different from forms 1 through 3. 

A universal principal-branch formula for transforming a nested power to an unnested 
power is 

{wf'y ^ {-ly w^\ (14) 

where 

Jo, ifarg(0) = 0, 

''"•=■( 7 (arg (u;^) — /3 arg (w)) (15) 

— ^^ ^^ -^ otherwise, 

V TX 

with short-circuit evaluation so that the "otherwise" result expression is not evaluated 
when the "if" test is true. 
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The transformation given by formulas (fT^ and ( TTSj) can be derived from the identities 

\p\ = (-l)--s(p)/-p ioip^O, (16) 

Iq^'f = |gr^ (17) 

Notice that the unit-polar factor (—1)^ is unit magnitude because arg(. . .) is always 
real, as are the rational numbers 7 and /3. Moreover, (—1)^ is piecewise constant with 
pie-shaped pieces emanating from w = because arg (^w^^ and f3 arg (w) have the same 
derivative with respect to w everywhere they are both continuous, and each of them has 
a finite number of discontinuities. 

An imaginary exponential e*'^'^ is an alternative to (—1)^, but it has the candidness 
disadvantage of introducing i into a factor that can be real and always is for the common 
case where the outer exponent 7 is a half-integer. 

If arg(O) is defined as 0, as it is in Mathematica, Maple, and Maxima, then we can 
define r more concisely and unconditionally as 

7(arg(u;^) -/3arg(ii;)) 
r := — ^ ^ -. (18) 

TV 

Proposition 2. Ifw> 0, then r = 0. 

Proof. When w = 0, r = follows immediately from expression (ITSl) . and 

u; > ^ arg (ty^) = A arg (w) = 0^7 (arg (ty^) — /3 arg (w)) /tt = ^ r = 0. D 

Proposition 3. If -1 < (3 < 1, then r = 0. 

Proof —1</3<1 =^ arg (w^) = (3 arg (w) =^ 7 (arg (w^) — (3 arg (w)) /vr = =^ t = 
0. D 

Proposition 4. 1/7 is integer, then (—1)^ = 1- 

Proof, arg (w^) is (3 arg (w) plus an even integer multiple of 2-7r. Thus when 7 is an 
integer, then 7 (arg (w^) — (3 arg (w)) /vr is an even integer, making r be an even integer, 
making (—1)^ = 1. D 



The simplification afforded by these three propositions should have already been ex- 
ploited with bottom-up default simplification, in which case [w^^ will have already been 
simplified to w^"' . If it isn't, then that is another opportunity to improve the system for 
very little effort J^i Thus, because {3 and 7 are explicit no n- zero rational numbers, without 
loss of generality this article assumes that w isn't known to be nonnegative, and that 
/3<— lor/3>l, and that 7 is non- integer. 

Using transformation (IT^ on every {w^^^ '' in W defined by ([1]) then collecting powers 
of —1 gives 

where cr is a simplified sum of terms of the form (fT5|) or flTSj) . 

The factor (—1)*^ is also unit magnitude with pie-shaped piecewise constant pieces 
because it is the product of such factors. This form has two great advantages over the 
other three forms: 



^°Do your computer algebra system's default and optional transformations de-nest (w") for such (3, 
7, and w declared non- negative? 
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• All of the exponents have been combined into a single unnested exponent. 

• Cancelable singularities are always completely canceled. 

Unfortunately this conies at the expense of a form that is usually bulkier than the other 
forms 

Simplification of individual piecewise expressions and combinations of such expres- 
sions is currently rather weak in most systems, but Carette [2] describes a canonical form 
for such expressions, so we can hope for improvement. In our case the piecewise expres- 
sions all have the same tests. Therefore we can add all of the Os together and add all of 
the expressions involving arg(. . .) together into a single piecewise function. For example. 



.2n3/2 / 3n4/3 




if aTgz = 0,\ I 
otherwise 



|(arg(z3) 



-3 arg z 1 



if arg2; = 0,' 
, otherwise 



/ 



-)■ 






^|2+|3-6 



I arg(^^^j + | arg(2^)-7arg2 



-1^ 



v 



if arg 2; = 0, \ 
otherwise 

/ 



z. (20) 



If arg(O) — )■ 0, then simplification of piecewise expressions isn't an issue here and 
the resulting exponent of —1 is simply (| arg {z"^) + | arg{z^) — 7 arg z) /vr. However, the 
result is not canonical either way, because starting with the equivalent canonical form 2, 



z^ (z^) 



3a1/3 



— > 



/ 



\ 



0, 

i arg(z2) + i avg{z^)-2a,Tgz 



if arg z = 0, \ 



otherwise 



(21) 



/ 



which has smaller magnitude coefficients. Thus for canonicality we could precede this 
transformation with a transformation to form 2. Equivalently we can adjust the coef- 
ficients of the arg (w^*^) and arg(w) analogous to how we adjusted exponents to arrive 
at form 2. This is preferable because it also canonicalizes expressions of form 4 that are 
entered directly or generated by the system. 

With pie-shaped pieces, ( — 1)°^ can always be expressed in the more candid canonical 



form 



Ci, if — 7r<argif; :: 9i^ 
C2, if 9i :: argw :: 62, 



v^Cm, otherwise. 
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(y'Y' 


-1 


r 
1 


if 5R (w;) > V 5R (w;) > A 53 (w;) > 


W 


otherwise; 






'1 


-5RH <53H <^{w), 


I 4\l/4 


— )■ < 


— Z 
-1 


li -^{w)<^ iw) <^{w), 


W 


5RH <^{w) <-^{w), 






i 


otherwise. 



where Ci through c^ are unit-magnitude complex constants, 9i through 9m-i are real 
constants in (— 7r,7r), and each instance of"::" is either "<" or "<"|^ Moreover: 

1. When w is real, then the positive and negative real axes are each entirely within one 
pie slice, enabling us to simplify (—1)'^ to one unconditional constant or piecewise 
expression of the form 

Ci, if w :: 0, 

C2 otherwise, 

where "::" is one of the comparison operators ">", ">", =, "<", "<", or "7^". 

2. For half-integers or quarter-integer fractional powers, (—1)^ can be expressed as a 
piecewise expression depending on the real and imaginary parts of w rather than 
arg(u;). For example, 

(22) 



(23) 



Notice that the right side of result (122|) is the definition of the Maple csgn function. 

Without an abbreviation such as csgn(. . .), Most implementers will probably want to 
avoid form 4 as a default even when arg (0) — )■ 0, because (— l)'^ is likely to be rather 
complicated nonetheless: 

1. It will probably contain complicated square roots and arctangents if the real and 
imaginary parts of w are given as exact numbers. 

2. It will probably also contain piecewise sign tests if given real and imaginary parts 
that are non-numeric, such as for w = x + iy with non-numeric real indeterminates 
X and y. 

3. It will probably contain radicals nested at least one deep if arg (w) is a simple 
enough rational multiple of vr. 

4. Otherwise it will contain perhaps bulky sub-expressions arg (w) and arg (^w^) - or, 
worse yet, expressions involving square roots, arctangents, piecewise sign tests, and 
sub-expressions of the form 3f? (w) and 9= (w). 

As espoused by Cor less and Jeffrey |1|, expression r can alternatively be defined in terms 
of the unwinding function k as: 

T:=2-fK{l3lnw). (24) 



-^^In a degenerate case, one or more of the pieces of pie might be a ray - very dietetic. 
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This is more concise than definition (IT5|) . but a function that computes unwinding num- 
bers isn't currently available externally in most computer algebra systems. Also, unless 
the system automatically transforms InO to — oo, as is done in Mathematica and Derive, 
then definition flM|) has the same disadvantages as using arg (. . .)o 

5.5 Simplifying mixtures of form 4 with form 1, 2 or 3 

If an expression contains a mixture of forms, then we should unify the forms to facilitate 
collection and cancelation. For example with arg(O) — )■ 0, the three expressions 

^^L, (25) 

z 

(_l)(a>-s(^')/2-arg^)A^ (26) 

1 if 3f? (w) > V 3fJ (ty) > A S (w) > 0, 
— 1 otherwise 

are equivalent. Therefore the result of any linear combination of them should transform 
either to or a multiple of one of them. The rewrite rules in the Appendix don't address 
this issue. 

In general it is easy to transform form ( l25l) to form (l26ll . which is only slightly more 
difficult to transform to either form (l25l) or form (1271) . 

6 Unimplemented extensions 

6.1 More semantic pattern matching for w 

The Mathematica pattern matcher is mostly syntactic rather than semantic, and the 
rules in the Appendix do almost no transformation of the radicand expressions w or any 
cofactors thereof. Thus recognition of opportunities relies mostly on the default trans- 
formations together with any optional transformations done by the user. Consequently, 
opportunities for the rules to simplify nested power products might not be recognized 
for radicands that aren't indeterminates. The rules work for most functional forms that 
have syntactically identical forms for the different instances of w, such as 

iLog[x^ {x + y)] j 2/3 

Log[xHx + y)] ' ^"^ [^' ^' + ^)] V-'^ [^' ^' + ^)]J ■ 

However the rules don't apply to all such functional form opportunities. For example. 



Cos[^] 



-)■ 



(Cosi^l^j'^/'Secf^] 



^^For TI-CAS, ln(0) — > undef. An error is inconveniently thrown by Maple for ln(0) and by Maxima 
for log(O). 
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because default simplification transforms Cos[^]^^ to Sec[^]. 

Even more opportunities are unrecognized when w is a sum. As an example of how 
to overcome this, the Appendix includes one extra rule that square-free factors radicands 
that are sums so that, for example. 



-> ^ 7^ y 



iz + i){iz + iyy 



z + l 

Factored over the integers or square-free factored form is a good choice for radicands 
for other reasons too, and these forms are canonical when the radicand is a rational 
expression. However, 

{z^ + 2Z+1) {{z + Iff" ^ (^ + 1)' ((^ + 1)')'^' ^ ((- + 1)')'^' , 

would require another rule that factors the cofactor of a power of a power of a sum. 
Then, perhaps we would want another rule to factor sums containing such radicands so 
that 

z^ {{z + 1)^)^/^ + 2z ((. + 1)^)^/^ + ((. + 1)^)^/^ -, {{z + 1)^)^/^ . 

It is impossible to implement equivalence recognition for all possible expressions w rep- 
resentable in general purpose systems, but it is worth expending a modest amount of 
execution time for default simplification and more time for optional transformations. 

The Appendix leaves most such opportunities unimplemented because the simplifi- 
cations described here are so fundamental and low level that they should be part of the 
built-in transformations. Good simplification of nested power products is more appro- 
priately built into a system rather than provided as an optionally loaded package that 
most users are unlikely to know about and load into every session. So rather than imple- 
menting a comprehensive package for one system, the intent of this article is to inspire 
implementers of all systems to improve some very fundamental transformations - at least 
to the extent that it can be done economically. 

6.2 Non numeric exponents 

Although not implemented in the rules of the Appendix, more generally the exponents 
for forms 1 through 4 can be Gaussian fractions or even symbolic, in which case we can 
still apply these transformations to the rational numeric parts of the exponents. For 
example, 

As another example, if a user has declared the variable n to be integer, then 

To some extent, the methods can also be extended to handle floating-point and sym- 
bolic real expressions for exponents a and P^. For example. 



w^-'"' 



(^1-234)3/2 ^ ^5.555y^;;T:^^ 



^2-.. .N3/2 ^ ^2 , 
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7 Summary 

This article: 

1. shows that many widely-used computer algebra systems have significant room for 
improvement at simplifying sub-expressions of the form w" (w^^) ^ ■ ■ ■ (w''") "; 

2. defines four different simplified forms with good properties; 

3. explains how to compute these forms; 

4. includes a demonstration implementation of form 3 via Mathematica rewrite rules. 
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Appendix: Mathematica rewrite rules for w" (w^^Y^ • • • (w^'^Y'' 

(* EXTRA SIMPLIFICATION DONE BEFORE ORDINARY EVALUATION: *) 

PreProductOfPowersOf Powers [(w_Plus)~(g_Rational /; ! IntegerQ [g] )] := 
Block [{squareFree = FactorSquareFree [w] } , 

squareFree"g /; Head [squareFree] =!= Plus]; 

PreProductOfPowersOfPowers [(w_"b_)"(g_Rational /; g <= -1 I I g >= 1)] : = 
w^ClntegerPart [g] *b) * (w"b)"FractionalPart [g] ; 

PreProductOfPowersOfPowers [(w_~b_)~(g_Rational /; g<=-l II g>=l) * w_~a_. * u_] := 
PreProductOfPowersOfPowers [w"(a+IntegerPart [g] *b) * (w"b)"FractionalPart [g] * u] ; 

PreProductOfPowersOfPowers [(w_"b_)"g_ * w_"a_. * u_. /; Sign [a] != Sign[b*g] && 
(Sign [a+b*Sign[g]] == Sign [b*(g-Sign[g] )] I I 

Min [Abs[a], Abs[b*g]] > Min [Abs [a+b*Sign[g]] , Abs [b*(g-Sign[g] )]] I I 
g == -1/2 && Min [Abs[a], Abs[b/2]] == Min [Abs [a-b] , Abs[b/2]])] : = 
w~(a+b*Sign[g]) * (w"b)"(g-Sign[g] ) * u; 

PreProductOfPowersOfPowers [f_[args ]] := 

Apply [f, Map [PreProductOfPowersOfPowers, {args}]]; 

PreProductOfPowersOfPowers [anythingElse_] := anythingElse ; 

(* EXTRA SIMPLIFICATION DURING ORDINARY EVALUATION: *) 

Unprotect [Times] ; 

(w_"b_)"g_ * w_"a_. * u_. /; Sign [a] != Sign[b*g] && 
(Sign [a+b*Sign[g]] == Sign [b*(g-Sign[g] )] I I 

Min [Abs[a], Abs[b*g]] > Min [Abs [a+b*Sign[g]] , Abs [b*(g-Sign[g] )]] I I 
Min [Abs [a], Abs[b*g]] == Min [Abs [a+b*Sign[g]] , Abs [b*(g-Sign[g] )]] && 

Abs[g] > Abs [g-Sign[g]] I I 
g == -1/2 && Min [Abs [a], Abs[b/2]] == Min [Abs [a-b] , Abs[b/2]]) := 
w"(a+b*Sign[g]) * (w~b)~(g-Sign[g] ) * u; 

(w_"bl_)"gl_ * (w_~b2_)~g2_ * u_. /; Sign[bl*gl] != Sign[b2*g2] && 

Abs[b2] > Abs[bl] && Sign [b2*(g2-Sign[g2] )] == Sign [bl*gl + b2*Sign[g2]] := 
(w"b2)"(g2-Sign[g2]) * ((w"(b2*Sign[g2] ) * (w-~bl)-gl) * u) ; 

Protect [Times] ; 

(* EXTRA SIMPLIFICATION DONE AFTER ORDINARY EVALUATION: *) 
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PostProductOfPowersOf Powers [w.^a. * (w_~b_)'~g_ * u_. /; IntegerQ [a/b]] : = 
PostProductOfPowersOf Powers [(w"b)"(g+a/b) * u] ; 

PostProductOfPowersOfPowers [w_"a_. *(w_"b_)"(g_Rational /; g<=-l II g>=l)*u_. 
/; llntegerQ [(a + b*IntegerPart [g] )/b]] := 
PostProductOfPowersOfPowers [ (u*w" (a+b*IntegerPart [b] ) ) * (w^b) ~FractionalPart [g] ] ; 

PostProductOfPowersOfPowers [f_[args ]] := 

Apply [f, Map [PostProductOfPowersOfPowers, {args}]]; 

PostProductOfPowersOfPowers [anythingElse_] := anythingElse; 

$Post = PostProductOfPowersOfPowers; $Pre = PreProductOfPowersOf Powers; 

Tables [I through [HI 



Table 1: Unflawed results of Appendix rewrite rules for 1st row x 1st column. 
Compare with Tables |2] through [7] 



|x 


(^r'" 


(.=)-'» 


/ 2\-l/3 


(wr" 


(w^f 


{w^f" 


w~^ 


1 


1 


1 


1 


w 


1 2\2/3 


«)7(«,2)l/3 


w5(«,2)l/3 


u;3(«,2)i/3 


u,(«;2)i/3 


(«,2)l/3 


W-' 


1 


1 


1 


1 


/ 2\2/3 


{w'f 


(«;2)i0/3 


(«,2)7/3 


(«,2)4/3 


(«;2)i/3 


W-' 


1 


1 


1 


w 


w iw^f" 


w' {w'f 


^5(^2)1/3 


«)3(^2)l/3 


«)(«)2)l/3 


1 9n1/3 


w° 


1 


1 


1 


i^^r 


(W^f 


(w^f 


iu,2)y^ 


(«,2)*/3 


(«,2)l/3 


w' 


1 


1 


w 


w {w'f' 


w' {w'f 


w' {w'f" 


^3(^2)1/3 


«,(«,2)l/3 


(w')"-' 


w' 


1 


1 


i.Y' 


{w-r' 


{wT' 


{wT'" 


(«,2)^/3 


(«,2)V3 


w^ 


1 


w 


»(»^)^" 


w' {w^f' 


W'' (^2)2/3 


w^w^f'"" 


U,(U,2)1/^ 


(«,2)V3 
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Table 2: Mathematica 8 default simplification for 1st row x 1st column, with flaw 



mbers 


Compare with Table [T 


1 
















i? 


{wT''' 


(«;2)-^/^ 


K)-"» 


i^r" 


(i/;^)^^ 


(l^2)«/^ 


w^^ 


1 


2 

6 


1 


2 

6 


1 


2 


1 2N2/3 


2 

3 

5 
6 


(«,2)^/3 
Ul3 


2 

3 

6 


r ) 


2 

3 

5 
6 


,^3(^2)7/3 


«,3(^2)4/3 


?«3(^2)l/3 


«,3 


Ul3 


w~'^ 


1 


2 

4 
5 
6 

7 


1 


2 

4 
5 
6 

7 


1 


2 

4 
6 


1 2N2/3 


2 

3 

4 
5 
6 


K) 
^2 


2 

3 
4 
5 
6 


/ 2N«/3 


2 

3 

4 
5 
6 


«,2(^2)7/3 


«,2(^2)4/3 


«,2(^2)l/3 


W; 


^2 


W-' 


1 


2 

5 
6 


1 


2 

5 
6 


1 


2 


/ 2x2/3 


2 

3 

6 


K) 


2 

3 

5 
6 


r ) 


2 

3 

5 
6 


1«(«)2)7/3 


«)(«)2)4/3 


i(;(«)2)i/3 


UI 


UI 


w^ 


1 


2 


1 


2 


1 


2 


/ 2n2/3 


2 


{w^f" 


2 


(^2)«/^ 


2 


(■!«2)''/3 


(«,2)4/3 


(«,2)l/3 


w' 


W 


2 

3 

5 
6 


w 


2 

3 
5 
6 


UI 


2 


1 2\2/3 

1/7 (ti; j 


2 


wiw'f 


2 

6 


1/; (i/;^)'/' 


2 

5 
6 


{w^yr-i 


(«,2)4/3 


(«,2)l/3 


w' 


1 


2 


1 


2 


(u.r'' 


2 


{w'f" 


2 


{wT" 


2 


iwY^' 


2 


(«;2)'l/3 


(t«2)l/3 


w^ 


«;3 


2 

3 
5 
6 

8 


«;3 


2 

3 
5 
6 

8 


U,3 


2 

3 
6 


W^ (|L;2)2/3 


2 


W' (^2)5/3 


2 

5 
6 


w^ {w^f^ 


2 

5 
6 


{w^f/^ 


(«,2)4/3 


(«,2)l/3 



Table 3: Mathematica 8 FullSimplify[.. 

bers. Compare with Table [T] 



,] for 1st row x 1st column, with flaw num- 



ix 


{^r"'' 


(^2)-4/^ 


(^2)-V3 


i^r' 


/ 2n5/3 


( 2n8/3 


w~' 


(u2)- 


1 
3 

5 


(.2)2/3 


1 
3 
5 


(.2)2/3 


1 
3 

5 


("T 


3 
5 


UI 




1 2n2/3 




u,s 


U,^ 


u5 


w^ 


(.2)1/3 


w~' 


(u2)^/^ 


1 
3 
4 
5 
6 


(.2)^/3 


1 
3 
4 
5 
6 


1 




1 




in^^r 




iw^r 




w^ 


U|6 


(u2)'l/3 


(U,2)V3 


w~' 


K) 


1 
3 

5 


(u,2)^^^ 

U>5 


1 
3 

6 


^ 2N2/3 

r ) 


1 
3 
5 
6 


w 




w {w'f 




w' iw'f' 




U.7 


u3 


(.2)1/3 


w' 


1 




1 




1 




(-^'^ 




{w'f 




{w'f 




(«;2)7/3 


(u;2)4/3 


(u2)l/3 


w^ 


u 


3 

5 


w 


3 

5 


w 




w (w^) 




w (w^f 


5 
6 


( 2n8/3 

w [w ) 


5 
6 


(«,2)7/3 


(u2)^/3 


(u2)V3 


w' 


1 




1 




/ 2n2/3 




(w^f 




( 2\8/3 




( 2n11/3 




(u,2)4/3 


(U,2)V3 


w'^ 


U,3 


3 

5 
6 
8 


w 




wiw'f 




w^ {wY' 




w" {w'f'- 


5 
6 


w' (w'f 


5 
6 


(«,2)7/3 


(U2)V3 
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Table 4: Derive 6 default simplify for 1st row x 1st column, with flaw numbers. 
Compare with Table [H 



i? 


{u,r''-^ 


(^.^)-/^ 


(-^)-^^^ 


/ 2\2/3 


(»^)"» 


(w^f 


w~^ 


«,9 


1 

3 


«;7 


1 

3 


K)- 


1 

3 

5 


r ) 


3 


(..)- 


3 

5 


W (W7^) 




«,5 


UI'^ 


U) 


W'^ 


1 




1 




1 




1 




/ 2\2/3 




{w'f 




(«,2)l°/3 


(«,2)^/^ 


(w^f^' 


(«,^)V3 


W-' 


(«,2)^/^ 


1 
3 


(«,2)^/^ 


1 
3 

5 


(»2)^/- 


1 
3 
5 


r ) 

UI 


3 

5 


w{w^f' 




w^ {w^f^ 




w"^ 


ui5 


«,a 


^0 


1 




1 




1 




(-^)'"' 




iw'f 




{w'f 




(«,2)^/3 


(^^^j-l/a 


(«,2)l/=' 


w^ 


/ 2N2/3 


3 

5 


/ 2N2/3 


3 

5 


/ „n2/3 


3 

5 


^ 2\2/3 




w^ {w^f^ 




w' (l/;2)2/3 




U.5 


«)3 


UI 


2 
W 


1 




1 




(-=)'^" 




(i/;^)^/^ 




{w'f 




/ 2\ll/3 




(«,2)4/3 


(«,2)Va 


W^ 


r ) 


3 

5 


(«,2)^/- 


3 

5 


w (w^f 




M;3 (^i;2)2/3 




w' {w^f 




w'' K)'/3 




«>3 


UI 



Table 5: TI-CAS 
Compare with Tab! 


3.1 default simplify for 1st row x 1st column, with flaw numbers. 
eE 




ix 


{wr'" 


{^r^'-' 


(»')-"» 


/ 2n2/3 


(^2)^/^ 


(.^2)«/=^ 


w-^ 


1 




1 




1 


2 


(u;2) 


2 

3 

5 
6 


(u2)^^^ 


3 

5 
6 


w {w'f 




u>7(u;2)i/=* 


u5(u.2)l/3 


U,3(u2)l/3 


«,3 


w 


_9 
W 


1 




1 




1 




1 




( 2n2/3 




( 2n5/3 




(u2)^«/=* 


(U,2)^/^ 


(u2)4/3 


/ ?\l/3 


W^^ 


1 




1 




1 


2 


/ 2^2/3 


2 

3 

5 
6 


wiw'fl^ 




W^ (^2)2/3 




u>5{u;2)i/=' 


u3(u2)l/3 


U(u2)l''3 


w 


w° 


1 




1 




1 




i^-r" 




iu^r'' 




(«;2)«/^ 




(u2)^/3 


(u,2)*/3 


(U2)V3 


w' 


1 




1 




w 


2 


( 2^2/3 


2 


w^w'f' 




w;^(w;2)^/^ 




^3(^2)1/3 


U,(u,2)l/3 


(U2)V3 


w' 


1 




1 




/ 2n2/3 




( 2n5/3 




{w^f 




/ 2\ll/3 




(u2)4/3 


(u,2)l/3 


w^ 


1 




u 




u^ 


2 

3 
6 


m;^ (w^)'/' 


2 


w' {w^f" 




w'^iw-'f' 




U(u2)l''3 


(u,2)l/3 


(U2)V3 
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Table 6: Maple 15 and Maxima 5.24 default simplification for 1st row x 1st column, 
with flaw numbers. Compare with Table [H 



4,x 


( 2\-7/3 


/ 2N-4/3 


/ 2\-l/3 


(-^)'^' 


/ 2\5/3 


/ 2\«/3 


w-'^ 


1 


2 

5 
6 


1 


2 

5 
6 


1 


2 




2 

3 

5 
6 


(.2)^/^ 


2 

3 

5 
6 




2 

3 
5 
6 


^3(^2)7/3 


«>3{«)2 )'*/=' 


wi{-uflfl'^ 


«>3 


U'3 


U'3 


w-^ 


1 


2 

4 
5 
6 


1 


2 

4 
5 
6 


1 


2 

4 
6 


K) 


2 

3 

4 
5 
6 


K) 


2 

3 

4 
5 
6 


r ) 


2 

3 
4 
5 
6 


^2(^2)7/3 


^2(^2)4/3 


^2(^2)1/3 


«)2 


«;2 


Ul2 


w~^^ 


1 


2 

5 
6 


1 


2 

5 
6 


1 


2 


r ) 


2 

3 

5 
6 


r ) 


2 

3 

5 
6 


r ) 


2 

3 
5 
6 


«)(«)2)7/3 


«)(«)2)4/3 


u>(u;2)l/^ 


UI 


w 


UI 


w^ 


1 


2 


1 


2 


1 


2 


/ 2n2/3 


2 


{w'f 


2 


[w^fl' 


2 


(«,2)7/3 


(«;2)4/3 


(«,2)i/3 


w' 


w 


2 

3 

5 
6 


Ul 


2 

3 
5 
6 


Ul 


2 


t. (t.2)'/^^ 


2 


wiw'f 


2 

5 
6 


1 2n8/3 


2 

5 
6 


(«;2)^/3 


(«,2)^/3 


(«;2)i/3 


w^ 


t«2 


2 

3 
4 
5 
6 

8 


«,2 


2 

3 

4 
S 
6 
8 


^2 


2 

3 
4 
6 


w^ {w'f 


2 

4 
6 


w^ {w^f^ 


2 

4 
5 
6 


w^ {w^f" 


2 

4 
5 
6 


(«,2)^/3 


(«,2)^/3 


(^2)1/3 


w' 


W^ 


2 

3 

5 
6 

8 


«,3 


2 

3 
5 
6 

8 


«,3 


2 

3 
6 


w' {w'f 


2 


W^ (^2)5/3 


2 

5 


w' {w'f 


2 

5 
6 


(«;2)7/3 


(«;2)4/3 


(«;2)l/3 



Table 7: Maxima 5.24 fullratsimp(...) and Maple 15 simplify(...) for 1st row x 
1st column, with flaw numbers. Compare with Table [T] 



|x 


{wr'" 


(.=)-"» 


/ 2\-l/3 


/ 2\2/3 


{w'f 


{w^f 


w~-' 


1 




1 




1 




/ 2\2/3 


3 

5 
6 


1 2N2/3 

r ) 


3 

5 
6 


w {w'f 




«,7(^2)l/3 


«;5(«,2)l/3 


«,3(u2)l/=' 


u, 


_9 

w 


1 


4 


1 


4 


1 


4 


1 2N2/3 


3 

4 
5 
6 


{^^r 




W' (^2)2/3 


4 


«,6(«;2)i/3 


^«4(^„2)l/3 


u;2(„,2)i/3 


VJ 


w~^ 


1 




1 




1 




( 2^2/3 


3 

5 
6 


wiw'f 




w^ {w'f 




«>5(u,2)l/3 


u3(lo2)l/3 


1«(«)2)l/3 


w 


w° 


1 


4 


1 


4 


1 




iror- 




w' {w^f' 


4 


w' {w'^f' 


4 


«,4(^2)l/3 


«)2(l„2)l/3 


(u2)V3 


w' 


1 




1 




w 




w iw'f 




w^ {w'f^' 




w'' {w^f 




«,3(^2)l/3 


U,(u,2)l/3 


(u2)l/3 


w' 


1 


4 


1 




u2 


3 
4 
6 


w' {w'f" 


4 


w' («;2)2/3 


4 


w'iw'f' 


4 


«,2(^2)l/3 


(.,2)1/3 


(u2)l/3 


w^ 


1 




w 




IfJ'^ 


3 
6 


w' {w^f 




w'' {w^f 




W^W^f^' 




u(u2)l/3 


i^-fr^ 


(«;2)i/3 
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Table 8: Unflawed results of Appendix rewrite rules for 1st row x 1st column. 
Compare with Tables 19] through [151 



4,x 


iw-r''' 


(w-r"' 


iw-'r''' 


{n^-'f 


/ -9x3/2 


/ -9x5/2 


w~^ 


w 


\Ji"" 


v^ 


v^ 


v^ 


v^ 


\l^ 


w 


«r^ 


tU'"' 


ui'^ 


w-^ 


1 


1 


v^ 


(i)''^ 


(i)='^ 


l^) 


iw' 


yiii^ 


w~^ 


w'^ 


w 


v^^ 


v^ 


v^ 


v^ 


75" 


\f^ 


■U) 


«;3 


UI^ 


w^ 


1 


1 


1 


v^ 


(^r 


(^)^''^ 


/ 1 \V2 


(i.r 


v^ 


w' 


w^' 


w^ 


w 


v/^- 


v^ 


v^ 


v^ 


yi^ 


v^ 


«) 


-u)-^ 


w^ 


1 


1 


1 


1 


v^ 


/ 1 \3/2 


i^r 


/ 1 N5/2 


/ 1 N3/2 


v^ 


w^ 


w"^ 


w'' 


W^ 


li) 


v^^ 


v^ 


\f^ 


v^ 


\f^ 


v^ 


«) 



Table 9: Mathematica 8 default simplify for 1st row x 1st column, with flaw numbers. 
Compare with Table |H1 



ix 


( -2N-5/2 


(»-)-'^ 


(»-)-''^ 


(^-)^/^ 


{w~^f" 


iw^^f^' 


w~^ 


1 


2 

3 
5 
6 

8 


1 


2 

3 

5 
6 

8 


1 


2 

3 
6 

8 


v^ 


2 


i^r 


2 

5 
6 


i^r 


2 

5 
6 


I 1 N5/2 •> 


(^)"^-^ 


y^-^ 


Ul3 


Ul3 


Ul3 


w~' 


1 


2 


1 


2 


^^ 


2 


/ 1 x3/2 


2 


(^) 


2 


( 1 ^^/^ 


2 


/ 1 x3/2 


\lii^ 


W-' 


1 


2 

3 
5 
6 

8 


1 


! 
i 
) 
i 

i 


1 


2 

6 

8 


\l^ 


2 


(^)^^^ 


2 

5 
6 


i^r 


2 

5 
6 


/ 1 N5/2 


/ 1 N3/2 


y^"' 


w 


UI 


w 


w^ 


1 


2 


1 


2 


1 


2 


v^ 


2 


(i)"^ 


2 


ii^r 


2 


(^)^^^ 


/ 1 \3/2 


v^ 


w' 


tu 


2 

5 
6 


W 


2 

5 
6 


UI 


2 


\/s^ 


2 


(i)''^- 


2 

3 
5 
6 


i^r^- 


2 

3 

5 
6 


{i,r 


/ 1 n3/2 


v5- 


w^ 


2 


2 

4 
5 
6 


«;2 


2 

4 
5 
6 


«,2 


2 

4 
6 


V^-' 


2 

3 
4 
6 


(^r«' 


2 

3 

4 
5 
6 


i^r-' 


2 

3 

4 
5 
6 


(^)"^ 


/' 1 N3/2 


v^ 


w^ 


U,3 


2 

5 
6 


ifj'^ 


2 

5 
6 


«,3 


2 


v^-^ 


2 

3 
6 


(^r-' 


2 

3 
5 
6 


i^r-'' 


2 

3 

5 
6 


/ 1 N5/2 


/ 1 N3/2 


VI^ 
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Table 10: Mathematica 8 FullSimplify[...] for 1st row x 1st column, with flaw 



inibers 


. Compare with Table [8 


















J,x 


in,-')-"'' 


(^.-r'"^' 


(»-)-"^ 


(»-)"^ 


(--^)^^^ 


(^-2)^/^ 


w~^ 


V^-' 


3 
6 


v/ii"" 




v^ 




v^ 




(^)^^^ 


5 
6 


«,3 


5 
6 


w 


«)3 


w-i 


W-' 


1 




1 




v^ 




(i)''^ 




(^r 




(i)"^ 




(i^r 


-Ji 


W-' 


•J^'-' 


3 
6 


v/i""' 


3 
5 
6 


^/^- 




v^ 






5 


1' 1 N5/2 


5 


w 


w 


w^ 


1 




1 




1 




v^ 




(^)^'^ 




(^r 




{^r 


i^r 


v^ 


w^ 


w 


5 
6 


w 


6 


w 




v/^- 




{^r- 


6 


(^r- 


6 


(^r 


(i,r 


v^ 


w' 


sji"-' 


1 
3 

4 


sj^'-' 


1 
3 

4 


v'^'^-' 


1 

3 

4 


^^»^ 


3 
4 


v^ 




(^r 




w^ 


J^-" 


1 
3 


J^'"' 


1 
3 


y^-^ 


1 

3 


V^-' 


3 


\/S^ 




{^r-^ 


3 



Table 11: Derive 6 default simplify for 1st row x 1st column, with flaw numbers. 
Compare with Table |8l 



|x 


(»-)-'^ 




(--)-"^ 


[w-r'' 


(--)''^ 


(^-)^/^ 


w~^ 


-'^li 


3 


\J^^ 




v^ 




■j^ 




v;:^ 




v;:^ 




w 


w-^ 


«,5 


«;7 


W-' 


1 




1 




v^ 




i^r 




(i.)^'^ 




(^) 




(M"' 


^i. 


w-^ 


»v^ 


3 


"V^^ 


3 


V^^ 




v^ 




v^ 




v^ 




w 


u,-^ 


i(;^ 


w' 


1 




1 




1 




Y ^ 




(^) 




/ 1 \5/2 




{^r 


(i,r 


V^ 


w' 


»Vi 


3 


»v* 


3 


»v* 


3 


V^^ 




v^ 




v;:^ 




w 


w'i 


w' 


1 




1 




1 




1 




v^ 




(^r 




(i)"" 


(i)"" 


(i)"" 


-ji. 


w^ 


»v^ 


3 


»'\/* 


3 


-Vi> 


3 


""V^ 


3 


V^^ 




^/i. 




tu 
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Table 12: TI-CAS and Maple default simplification for 1st row x 1st column, with 
flaw numbers. Compare with Table [HI 



ix 


( -2N-5/2 


^ -2N-3/2 


iw-r''' 


/ -2n1/2 


^ -2n3/2 


/ -2n5/2 


w~^ 


1 


2 

3 

5 
6 
8 


1 


2 

3 

5 
6 

8 


1 


2 

3 

5 
6 
8 


v^ 


2 


( 1 VV2 


2 

5 


/ 1 \5/2 
tu3 


2 

5 




o/ 1 \3/2 


-^y^ 


Ul3 


w-^ 


1 


2 

3 
4 

5 
6 
8 


1 


2 

3 

4 
5 
6 
8 


1 


2 

3 

4 
5 
6 
8 


v;;^ 


2 

4 
6 


i^r 


2 

4 
5 
6 


i^r 


2 

4 
5 
6 


-\i,r 


o/ 1 \3/2 


V w^ 


«;2 


w^ 


w'^ 


W-' 


1 


2 

3 

5 
6 
8 


1 


2 

3 

5 
6 

8 


1 


2 

6 

8 


v^ 


2 


w 


2 

5 
6 


w 


2 

5 
6 


/ 1 N5/2 


,' 1 N3/2 


"-y^ 


w 


w^ 


1 


2 


1 


2 


1 


2 


y-^ 


2 


i^r 


2 


( 1 x5/2 


2 


/ 1 \5/2 


/ 1 ^1^/2 


v5 


w^ 


w 


2 

5 
6 


If; 


2 

5 
6 


UI 


2 


"v/^ 


2 


-{^r 


2 

3 

5 
6 


»(^)"^ 


2 

3 

5 
6 


( 1 \5/2 


/ 1 \3/2 


v^ 


w^ 


2 

If; 


2 

4 
5 
6 


«;2 


2 

4 
5 
6 


«,2 


2 

4 
6 


»v^ 


2 

3 

4 
6 


-= iM"'' 


2 

3 

4 
5 
6 


»M^)"' 


2 

3 

4 
5 
6 


I' 1 \5/2 


(l^)"^ 


Vi^ 


w^ 


«,3 


2 

5 
6 


UI^ 


2 

5 
6 


«,3 


2 


»v* 


2 

3 


«' (^r 


2 

3 

5 
6 


■-•Hi^r 


2 

3 

5 
6 


(^)^^^ 


/ 1 \3/2 


V^ 



Table 13: Maxima 5.24 default simplify for 1st 
Compare with Table [HI 



row X 1st column, with flaw numbers. 



i? 


{w-r"' 


(»-)-»'^ 


(»-)-'■' 


(--)"^ 


/ -2\3/2 

{w 2) 


(l.-^)^/^ 


w~^ 


K)^/2 


1 

2 

3 
5 
6 


r ) 


1 

2 

3 
5 
6 


1 


2 

3 
6 

8 


VliT^ 


2 


1 


1 
2 

5 
6 


1 


1 

2 

5 
6 


U'3 


U'3 


y^-^ 


«>3 


^3(^2)3/2 


t«3(«,2)5/2 


w-^ 


K) 


1 
2 

3 
4 

5 
6 


/ 2\3/2 

r ) 

Ul2 


1 
2 

3 

4 
5 
6 


1 


2 

3 

4 
6 
8 


v^ 


2 

4 
6 


1 


1 
2 

4 
5 
6 


1 


1 
2 

4 
5 
6 


Ul2 


v^-^ 


«)2 


t«2(«,2)3/2 


«)2(^2)5/2 


W-' 


/ 2\5/2 


1 
2 
3 

6 


/ 2\3/2 


1 
2 

3 

5 
6 


1 


2 

6 

8 


y^ 


2 


1 


1 
2 

5 
6 


1 


1 
2 

5 
6 


«) 


UI 


V^'^ 


w 


U>(u)2)3/2 


u;(u;2)5/2 


w' 


(u-^)"^ 


1 
2 


in^^r 


1 
2 


1 


2 


v^ 


2 


1 


1 
2 


1 


1 
2 


75" 


(«,2)3/2 


(u2)5/2 


w^ 


/ 2n5/2 


1 
2 

5 
6 


w {w'f 


1 
2 

5 
6 


ii) 


2 


V'i'«' 


2 


w 


1 
2 
3 
5 
6 


«) 


1 
2 
3 
5 
6 


v^ 


(«,2)3/2 


(«,2)5/2 


w' 


w' {w'f 


1 
2 

4 
5 
6 


w' {wY' 


1 
2 

4 
5 
6 


«,2 


2 

4 
6 


v/A-^ 


2 

4 
6 


«,2 


1 
2 

3 

4 
5 
6 


^2 


1 
2 

3 

4 
5 
6 


75" 


(«,2)3/2 


(«;2)5/2 


w^ 


w^ (w^f 


1 
2 

5 
6 


w' {wY' 


1 
2 

6 


«,3 


2 


V'i'-' 


2 

3 

6 


«,3 


1 
2 

3 

5 
6 


w^ 


1 
2 

3 

5 
6 

8 


v^ 


(«,2)3/2 


(w^f^ 
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Table 14: Maxima 5.24 fullratsimp(...) for 1st row x 1st column, with flaw numbers. 
Compare with Table |H1 



4,x 


{w~^r"' 




(«;-)-/^ 


iw-^f' 


(^-)^/^ 


{w ') 


w~^ 


wyw"^ 


1 

6 


W 


1 

2 

6 


1 


2 

3 
6 

8 


v^ 


2 


1 


1 

2 

6 


1 


1 

2 

6 


y^-^ 


«,3 


w^^/w'^ 


ui'^Vui^ 


W-' 


W^\/%lP' 


1 
4 
6 


\/w^ 


1 
2 
6 


1 


2 

3 
4 
6 
8 


v^ 


2 

4 
6 


1 


1 
2 

4 
6 


1 


1 
2 

4 
6 


V^-^ 


li? 


t/J'^Vt/J^ 


ui^V^/j^ 


w~' 


yj'i^yji 


1 
6 


wy/w'^ 


1 
2 
6 


1 


2 

6 

8 


v^ 


2 


1 


1 
2 
6 


1 


1 
2 

6 


y^"' 


-U) 


w^^/w'^ 


t/j^V^^ 


w^ 


w'^Vw'^ 


1 
4 


iiP'y/vfi 


1 
2 

4 
6 


1 


2 


v^ 


2 


1 


1 
2 

4 
6 


1 


1 
2 

4 
6 


v^ 


UI^V^"^ 


UI^V^"^ 


w' 


w^y/w"^ 


1 
6 


w^y/w'^ 


1 
2 
6 


«) 


2 


\[i^ 


2 


1 


1 
2 

6 


1 


1 
2 
6 


v^ 


w^Jw'^ 


w^yjw'^ 


w^ 


yj^y/yj'i 


1 
4 
6 


w'^y/w'^ 


1 
2 

4 
6 


UI^ 


2 

4 
6 


\/>^ 


2 

4 
6 


1 


1 
2 

4 
6 


1 


1 
2 

4 
6 


v^ 


ui^Vui^ 


w^ 


nP yjiip' 


1 
6 


w^y/w'^ 


1 
2 
6 


w^ 


2 


v^-^ 


2 

3 
6 




1 
2 

6 


1 


1 
2 

6 


v^ 


ui-\/«)2 



Table 15: Unflawed results of Maple simplify(...) for 1st row x 1st column - a 
variant of form 4. 
Compare with Table [HI 



ix 


( -2N-5/2 


( ~2N-3/2 


(^-2)-V2 


(--)"^ 


{^^'f 


{n^^-f^ 


w~^ 


w^'csgn (i) 


csgn {£) 


csgn(i;) 
«)2 


«)4 


«;6 


csgn(i;) 


t«8 


w~^ 


w^csgn (i) 


u;csgn(i) 


UI 




csgn(i^) 
«)5 


„,7 


W-' 


ui-^csgn (i) 


w^'csgn (^) 


csgn (i) 


«)2 


csgn(i^) 


csgn(i;) 


t«6 


w^ 


w^csgn (^) 


w^'csgn (^) 


wcsgn(^) 


U) 


csgn(i^) 


-gn(^) 


«)3 


,«5 


w' 


w^'csgn (i) 


u^^csgn (^) 


w^csgn (i) 


csgn (i) 


csgn(i) 
iii2 


csgn(i;) 


W 


w^csgn (i) 


w^csgn (^) 


w^'csgn (i) 


wcsgn(i) 


csgn(^) 


-gn(^) 


U) 


„,3 


W^ 


^'csgn (i) 


w^'csgn (^) 


ty'^csgn (i) 


ty^csgn (^) 


csgn (i) 


csgn(^) 
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Table 16: Simplification of w /yw'^, with flaw numbers 



system 


transformations 


input ^ 




Appendix rewrite rules 


11) 




Mathematica 


default 


w 


2,6,8 


Mathematica 


FullSimplify(...) 


w 


6, 8 


Derive 


default 


111 




TI-CAS 


default 


w 


2,6,8 


Maple 


default 


w 

/ 9 


2,6,8 


Maxima 


default 


w 


2,6,8 


Maxima 


fullratsimp(...) 


ID 
/ 9 


6, 8 


Maple 


simplify (...) 


csgn (w) 





Table 17: Simplification of v 


WyW - (_l)(V2)(arg(«,2)_2arg{«;))A^ ^^^^ ^^^ numbcrS. 




system and transformation 


result and flaw numbers 


Appendix rewrite rules 


^(-2ArgH+Arg w^ )/7t , y^ 


2 


Mathematica default 


•(-2ArgH+Arg w'^ )/-k , y^ 
' w 


2 


Mathematica FullSimplify(...) 


•(-2ArgH+Arg w^ ')/tt . Vw'^ 
' w 


2 


Derive default 


"^ IPfty-O ( l)(PHASE(«;2)-2PHASE{«,))/{2^)\ 


2 


TI-CAS, default 


0, w = 

^ '"*'] (anglc(«)2)„2anglc(«;))l/2 ^ 

V^ g ^ else 

w 


2 


Maple default 


argument 1 UJ j — 2 arguinent(z) 


2 


Maxima default 


atan2{sin(2 carg(Tii)),cos(2 carg(ii.'))) — 2 carg(ii.') 


2 


-f^-(-^) 


Maxima fullratsimp(...) 


carg(i(;) atan2(sin(2 carg(i[f) ),cos(2 carg(i^))) 


2 


yup-ly — l) T _^(_ij ^TT 


. . <:arg(ti.) 


W{-V) T 


Maple simplify (...) 


argument 1 vj 1 —2 arguincnt(2) 

csgn {vS) — (— 1)^ 


2 
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